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ABSTRACT 


Direct  application  of  two-phase  heat  transfer  in  the  liquid  cooling  of  electronic 
components  in  flouorinated  hydrocarbons  (FC-72),  is  severely  inhibited  by  the  excessive 
amount  of  superheat  required  to  initiate  nucleate  boiling.  To  conduct  an  experimental 
study  of  nucleate  pool  boiling  of  FC-72,  an  experimental  test  chamber  was  constructed. 
This  chamber  utilized  five  horizontal  platinum  wires  of  0.05  mm  diameter  spaced  2.0 
cm  vertically  from  each  other.  The  lowest  wire  was  progressively  heated  from  the  na¬ 
tural  convection  region  through  nucleate  boiling,  and  a  study  was  made  on  the  effects 
of  the  boiling  wake  plume  on  the  heat  transfer  rate  of  the  upper  wires. 
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I.  IN  FRODUCTION 


A.  HISTORICAL 

Ever  since  the  devclopmciit  of  the  first  practical  solid  state  integrated  circuit  in  1959, 
there  has  been  a  continuous  emphasis  on  packaging  an  ever  higher  density  of  logic  cir¬ 
cuits  on  a  single  substrate.  Corresponding  with  this  increase  packaging  density  is  an 
increasing  power  dissipation  requirement  along  with  an  increase  difficulty  in  the  thermal 
management  capabilities.  Initially,  thermal  management  techniques  consisted  of  noth¬ 
ing  more  than  a  heat  sink  and  natural  connection  air  cooling.  This  naturally  progressed 
into  heat  sinks  with  forced  convection  air  cooling,  and  with  the  construction  of  large 
mainframe  computers  and  first  generation  super-computers  the  use  of  indirect  liquid 
cooling  was  initiated.  The  next  practical  method  for  heat  removal  is  direct  liquid  cooling 
or  immersion  cooling  of  microelectronic  components  using  dielectric  liquids. 

Although  direct  liquid  cooling  has  been  successfully  employed  for  nearly  40  years  in 
large  hi-powered  microwave  equipment,  klystron  tubes,  travelling-wave  guide  tubes,  and 
power  transistors  (Ref  1  :  p.l2],  only  recently  has  it  been  employed  for  use  in  tne  ther¬ 
mal  management  of  microelectronics.  The  first  comp:’tcr  systems  to  successfully  employ 
immersion  cooling  of  microelectronics  were  the  CR.'\Y-2  and  the  ETA-10  supercom¬ 
puter  (Ref  2  :  p.  8).  Today  the  research  in  direct  immersion  cooling  is  continuing  with 
particular  emphasis  on  utilizing  dielectric  liquids  as  a  cooling  medium,  in  particular  the 
prcfluorinated  FC  liquid  family  of  the  3M  corporation  [Ref  3). 

B.  BACKGROUND 

The  advantages  offered  by  direct  immersion  cooling  is  self  evident  from  Figure  1 
By  utilizing  forced  convection  direct  liquid  cooling  with  FC-77,  the  CR/\V-2  has  been 
able  to  reduce  it's  required  volume  by  an  order  of  magnitude  as  compared  to  it's  indirect 
liquid  cooled  competitors.  This  was  being  achieved  in  the  CR.‘\Y-2  while  at  the  same 
time  maintaining  the  same  power  dissipation  of  the  competition.  (Ref  1]. 
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Figure  1.  Power  consumption  vs  Volume:  Of  present  computers 

As  previously  stated,  the  CRAY-2's  packaging  density  was  achievable  due  to  direct 
liquid  immersion  cooling,  however,  it  did  use  a  forced  convection  cooling  scheme.  A 
thermal  management  scheme  that  relies  on  boiling  heat  transfer  would  provide  an  even 
more  efficient  heat  transfer  capability,  and  has  been  the  attention  of  considerable  re¬ 
search  in  recent  years.  Unfortunately,  these  prefluorinated  dielectric  fluids  possesses 
anomalies  which  have  held  back  their  widespread  use  and  acceptance  of  direct  immersion 
cooling. 

In  comparison  to  water  and  many  other  coolants,  these  prefluorinated  dielectric 
liquids  posses  a  lower  value  for  their  thermal  conductivity,  and  heat  of  vaporization, 
they  also  have  an  extremely  small  surface  tension  which  translates  into  a  near  zero  wet¬ 
ting  angle  and  therefore,  a  high  wettability  on  most  surfaces.  This  last  feature  results  in 
the  elimination  of  many  potential  nucleation  sites  which  in  turn  results  in  a  incipience 
super  heat  excursion  which  results  in  a  very  high  wall  temperature. 

In  this  study,  heat  transfer  enhancement,  and  therefore  a  surface  temperature  re¬ 
duction,  due  to  bubble  pumping  in  FC-72  will  be  investigated. 

C.  PREVIOUS  WORK 

An  extensive  amount  of  research  has  been  conducted  in  the  boiling  regime  of 
dielectric  liquids.  The  most  helpful  of  the  literature  used  will  be  noted. 
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You,  Simon  and  Bar-Cohen  [Ref.  4  and  Ref  5  ]  conducted  an  analyses  on  the  ef¬ 
fects  of  pressure,  subcooling,  and  dissolved  gas  content  on  the  boiling  heat  transfer,  and 
the  boiling  incipience  of  dielectric  fluids.  In  their  experiment  they  used  FC-72  as  the 
working  fluid  and  a  O.l^m  thin  film  of  platinum  on  a  0.51mm  diameter  quartz  cylinder, 
with  the  platinum  film  acting  as  both  heater  and  temperature  sensor.  They  concluded 
that: 

1.  An  increase  in  pressure  of  56kPa  reduces  wall  superheat  by  approximately  S'C  in 
the  fully  developed  nucleate  boiling  region. 

2.  Subcooling  has  little  effect  on  the  fully  developed  nucleate  boiling  region. 

3.  Nucleate  boiling  hysteresis  is  observed  for  both  the  subcooled  and  gas-saturated 
cases.  The  hysteresis  decreases  as  the  gas  content  increases. 

4.  Incipience  bubble  radius  computed  for  FC-72  was  0.05-0.  Ijum. 

5.  Effects  of  dissolved  gases  in  the  fluid  is  seen  at  low  content  values  as  an  increase 
in  run  to  run  variability. 

6.  The  effect  of  subcooling  on  pool  boiling  incipience  is  apparently  small. 

7.  A  wide  variation  in  incipience  superheat  values  for  nominally  identical  runs  were 
observed. 

You,  Simon,  Bar-Cohen  and  Tong  [Ref  6]  conducted  an  experimental  investigation 
of  boiling  in  R-113.  The  experiment  was  conducted  using  0.13mm  diameter  Chromel 
wires  and  0.5 1  mm  diameter  platinum  thin  film  heater.  The  following  was  found. 

1.  A  Scanning  Electron  Microscope  (SEM)  analysis  was  conducted  on  the  surface  of 
the  thin  film  heater,  in  which  the  chromel  wire  surface  was  found  to  contain  a 
higher  density  of  nucleation  sites  compared  to  the  thin  film  platinum  heater. 

2.  The  platinum  heater  surface  was  found  to  have  a  higher  incipience  superheat  than 
the  chromel  wires ,  73°  C  versus  38°  C  maximum.  This  tends  to  correlate  well  with 
item  1. 

3.  Again,  a  wide  variation  in  incipience  superheat  values  for  nominally  identical  runs 
were  experience. 

Marto  and  Lepere  [Ref  7]  conducted  experiments  with  surface  enhancements  for 
improving  pool  boiling  performance.  Marto  and  Lepere  using  FC-72  and  R-113  as 
working  fluids  compared  the  performance  of  three  commercially  available  boiling  en¬ 
hancement  surfaces  to  plain  copper  tubing.  In  each  case  the  enhanced  surface  resulted 
in  a  lower  wall  surface  temperature  compared  to  the  plain  copper  tube.  However,  the 
required  incipience  superheat  for  nucleation  remained  unaltered. 

Bergles  and  Kim  [Ref  2]  analyzed  methods  to  reduce  the  temperature  overshoot  of 
immersion  cooling.  Their  method  used  an  active  bubble  generator  in  the  fluid  located 
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below  the  simulated  micro  chip.  They  used  the  bubble  generator  to  activate  dormant 
and  extinct  nucleation  sites.  They  obtained  a  reduction  of  incipience  superheat  required 
from  22-33°  C  to  8-12°  C  in  a  working  fluid  of  R-113. 

Maddox  and  Mudawar  [Ref.  8]  studied  the  effect  of  subcooling  and  surface  aug¬ 
mentation  on  values  of  Critical  Heat  Flux  (CHF).  They  found  that  the  CHF  of  FC-72 
can  achieve  a  level  greater  then  360  Watts/cm^  Avith  combined  use  of  subcooling  and 
surface  augmentation. 

Danielson,  Tousignat  and  Bar-Cohen  [Ref.  9]  conducted  an  analysis  of  a  large 
number  of  commercially  available  prefluorinated  inert  fluids.  The  method  of  investi¬ 
gation  used  a  0.25nim  diameter  platinum  wire  heating  element,  with  the  working  fluid 
at  saturation  conditions.  Important  points  reported  are: 

1.  The  critical  heat  flux  of  FC-72  was  reported  as  approximately  20.3  Watts/cm^ 

2.  The  nucleate  boiling  characteristic  of  highly  wetting  fluids  can  be  predicted  using 
correlations  developed  for  water  and  other  fluids,  but  are  highly  dependent  of  the 
heating  surface  condition. 

3.  Single  phase  convection  data  are  predicted  accurately  by  a  correlation  developed 
by  Kuehn  and  Goldstein.[Ref.  10] 

4.  The  superheat  excursion  at  incipience  boiling  always  appears,  but  the  magnitude 
varies  widely. 

5.  The  apparent  radii  of  nucleation  bubbles  were  calculated  to  be  0.1  to  0.2/tm. 

Kuehn  and  Goldstein  [Ref.  10]  has  obtained  a  correlation  of  natural  convective  heat 
transfer  for  infinite  cylinders.  This  correlation  was  used  in  this  thesis  data  analysis  due 
to  the  low  rayleigh  numbers  obtained  in  the  experimental  runs. 

From  the  3M  Product  Manual  [Ref  3]  it  can  be  found  that  the  highly  wetting 
dielectric  fluids  have  a  high  solubility  for  gases.  FC-72  can  contain  up  to  48%  air  by 
volume  compared  to  1.9%  air  for  water.  This  is  the  necessity  of  thoughly  degassing  the 
fluid  prior  to  the  experimental  runs. 

D.  OBJECTIVES 

The  objective  of  this  work  was  to  investigate  the  effects  of  the  boiling  wake  from  a 
heated  wire  upon  the  heat  transfer  rate  from  wires  placed  within  the  wake.  This  work 
was  carried  out  using  the  dielectric  fluid  FC-72. 
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11.  EXPERIMENTAL  APPARATUS 


A.  DESCRIPTION  OF  COMPONENTS 

The  experimental  apparatus  used  in  this  study  is  a  redesign  of  an  experimental  ap¬ 
paratus  used  in  a  previous  thesis.  Emphasize  on  this  design  was  to  minimize  systematic 
errors,  minimize  uncertainty,  and  to  correct  deficiencies  to  ensure  an  overall  increase  in 
the  data's  creditability  over  the  previous  apparatus,  [Ref  1 1]. 

In  Figure  2  an  overall  view  of  the  experimental  setup  is  shown,  identifying  all  major 
components  that  were  utilized.  Figure  3  is  a  photograph  of  the  actual  experimental 
chamber. 

The  apparatus  chamber  consists  of  a  box  constructed  of  0.75  inch  thick 
polycarbonate  sheets  with  a  0.25  inch  thick  aluminum  coverplate  fitted  with  a  rubber 
0-Ring  seal.  Inside  the  box  are  two  bulk  heaters  to  raise  the  fluid  to  saturation  tem¬ 
perature  for  degassing  purposes,  and  to  maintain  the  bulk  fluid  temperature  near  con¬ 
stant.  The  aluminum  cover  plate  is  mounted  with  three  thermo-electric  coolers,  and 
serves  as  a  condenser,  where  by  the  vapor  condenses  and  returns  to  the  chamber  by 
gravity.  An  insert  board  is  held  in  the  chamber  on  which  are  mounted  the  test  wires, 
tension  springs,  thermocouples,  and  voltage  pickofTs. 

1.  Test  Chamber 

The  interior  dimensions  of  the  chamber  are  6.0  inches  wide,  by  2.375  deep,  by 
6.25  inches  high  ,  as  shown  in  the  isometric  view  of  Figure  8.  The  entire  chamber  was 
assembled  using  IPN  #16  Acrylic  bonding  adhesive  and  a  bead  of  Di-Metheyl  Ethelene 
was  used  to  seal  all  seams  along  the  interior  surface  to  ensure  liquid  tightness. 

To  permit  access  for  instrumentation  and  power  leads,  and  as  a  means  to  ensure 
pressure  equalization,  a  0.75  inch  diameter  and  0,25  inch  diameter  hole  were  placed  into 
the  rear  face  of  the  chamber  at  the  lower  right  and  upper  right  comer  respectively,  as 
shown  in  Figure  4.  Mounted  securely  around  these  access  holes  are  lengths  of  tygon 
tubing  of  approximately  18  inches  in  length,  with  the  opened  ends  suspended  well  above 
the  chamber. 

On  the  interior  side  walls,  slots  are  milled  to  permit  the  placement  and  removal 
of  an  instrumented  wire  board  on  which  the  experiments  are  conducted.  Figure  6  is  a 
photograph  of  the  insert  board,  while  a  schematic  drawing  is  shown  in  Figure  7. 
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2.  Aluminum  Cover  Plate 

A  0.25  inch  thick  by  7.00  inch  long  by  3.25  inch  wide  aluminum  plate  was  uti¬ 
lized  for  use  as  a  cover  plate  for  the  chamber,  and  as  a  condenser  surface  for  the  FC-72 
vapor.  Mounted  on  the  interior  face  are  two  thermocouples  to  monitor  the  condenser 
surface  temperature,  as  shown  in  Figure  9.  On  the  exterior  face  of  the  plate,  three 
thermo  electric  coolers  were  installed  for  heat  removal.  To  ensure  a  proper  seal  between 
the  aluminum  plate  and  the  box  structure  ,  a  shallow  0-Ring  groove  was  machined  into 
the  interior  face  and  an  0-Ring  of  0.0624  inch  diameter  was  installed.  To  adequately 
secure  the  cover  plate  to  the  chamber  box,  eight  stainless  steel  screws  of  1/16  inch  di¬ 
ameter  are  used  to  compress  the  0-Ring  against  its  mating  surface. 

As  mentioned  previously,  three  solid  state  thermoelectric  cooling  devices  from 
the  MELCOR  corp  were  employed  to  cool  the  cover  plate  condensor.  They  w  ‘re  nom¬ 
inally  operated  about  the  0.4  amps  and  2.0  volts  operating  point,  from  [Ref  12]  their 
operation  curves  and  technical  data  is  available. 

3.  Heaters 

As  a  means  of  degassing  the  FC-72  liquid  and  for  bulk  tem.perature  control,  two 
strip  heaters  of  125mm  x  1 1  mm  were  employed.  These  strip  heaters  were  secured  to  the 
bottom  interior  surface  of  the  chamber  with  epoxy,  and  were  orientated  lengthwise  in  the 
chamber  such  that  one  heater  strip  will  be  in  front  of  and  one  heater  strip  behind  the 
insert  board.  To  raise  the  bulk  fluid  temperature  to  saturation  temperature  and  for 
degassing  purposed,  both  heaters  are  on  line  and  operated  at  1.8  amps  and  20  volts. 
Once  degassing  is  complete,  the  forward  strip  heater  is  electrically  isolated  and  the  rear 
heater  is  run  at  0.9  amps  and  20  volts.  This  will  maintain  the  forward  portion  of  ex¬ 
perimental  chamber  at  near  saturation  temperature  with  minimum  stratification  of  fluid, 
and  maintain  liquid  in  a  degassed  state  since  the  fluid  is  still  boiling  from  the  heater 
strip's  surface. 

4.  Thermocouples 

A  total  of  Six  copper-constantan  type  thermocouples  with  a  wire  diameter  of 
0.05  inches  (0.127  mm)  were  utilized  for  temperature  measurements.  Two  thermocouples 
were  bonded  to  the  underside  of  the  aluminum  condensor  plate  utilizing  OMEGA  Bond 
type  101  thermally  conductive  epoxy,  four  thermocouples  were  installed  on  the  insert 
board,  and  one  was  used  for  measuring  the  ambient  air  temperature. 

5.  Insert  Board 

An  insert  board  was  designed  to  mount  into  two  0.250  inch  wide  grooves  that 
were  milled  into  the  chamber  side  walls.  Onto  this  insert  board  are  mounted  four 
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thermocouples  to  measure  bulk  fluid  temperature,  five  platinum  wires  of  nominally  0.002 
inch  (0.05mm)  diameter  by  4.0  inch  long  for  conducting  the  actual  boiling  experiments. 
Fi''e  springs  are  used  to  maintain  the  platinum  wire  tension  so  no  sagging  develops  due 
to  thermal  expansion  of  the  wires.  In  Appendix  A  the  piatinum  wire  calibrations  and 
specifications  are  shown  in  detail. 

B.  INSTRUMENTATION 

In  the  design  and  assembly  of  this  experimental  apparatus,  five  channels  were  wired 
and  instrumentation  provided  to  measure  voltage  drops  at  select  pouits.  Each  of  these 
channels  had  a  2.0  ohm  precision  resistor  wired  in  series  with  ^  platinum  wire,  where 
the  voltage  drop  across  the  resistor  was  used  to  obtain  the  current  flow  through  the 
platinum  wire.  The  voltage  drop  measured  across  the  platinum  wire  is  used  to  obtain 
the  wire  resistance,  power  dissipated  and  the  heat  flux  of  the  platinum  wire.  The 
platinum  wire  temperature  is  obtained  from  its  resistance  vs  temperature  calibration 
curve  obtained  prior  to  experimental  runs.  Instrumentation  was  additionally  used  to 
obtain  the  bulk  fluid  temperature,  which  will  be  required  for  calculating  fluid  properties. 

1.  Platinum  Wire 

Each  of  the  five  platinum  wires  were  run  in  series  with  a  2  ohm  (±  1%)  precision 
resistor  of  25  watt  power  rating.  The  following  table  gives  the  platinum  wire  and  resistor 
combinations  which  were  used. 


Table  1.  MEASURED  RESISTANCE  VALUES  FOR  2  OHM  RESISTORS 


Wire  Chan¬ 
nel 

Measured  Resistance,  20  Sample 
average 

20  Sample  Std.  Dev. 

1 

1.9994a 

0.0009a 

2 

2.0076Q 

0.0005a 

3 

1.5998Q 

0.0003a 

4 

2.00 18Q 

0.001 3a 

5 

2.0025a 

0.0008a 

A  2.0  ohm  precision  resistor  was  chosen  over  the  previous  design's  0.1  ohms 
precision  resistor  in  order  to  reduce  the  calculated  currents  uncertainty,  and  there  by 
reduce  the  calculated  platinum  wire  temperature  uncertainty,  [Ref  11;  p.  25]. 
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2.  Power  Supplies 

The  following  power  supplies  were  utilized  and  their  primary  function  listed. 

1.  Hewlett  Packard  6214C  O-lOvolts/0-1  amps:  Used  to  power  the  thermal  electric 
coolers. 

2.  Hewlett  Packard  6214C  1-lOvolts  /  0-1  amps:  Used  to  power  the  constant  heat  flux 
platinum  wire. 

3.  WP711  0-40  volts  /  0-1  amps :  Used  to  power  the  var>’ing  heat  flux  platinum  wire. 

4.  Hewlett  Packard  6286A  0-20volts  /  0-10  amps  :  Used  to  power  Heater  Strips  for 
degassing  and  bulk  temperature  maintenance. 

5.  Hewlett  Packard  6289 A  0-40  volts  /  0-1.5  amps :  Used  to  power  an  auxiliary  illu¬ 
mination  system  or  an  augmenting  cooling  fan  for  the  thermal  electric  coolers. 

3.  Acquisition  Unit 

To  obtain  the  necessary  measurements,  a  Hewlett  Packard  3852A  Data  Acqui¬ 
sition  System  was  employed,  utilizing  a  Hewlett  Packard  300  series  desktop  computer. 
The  HP3852A  utilized  the  following  utility  modules  inserted  into  it's  chassis. 

1.  Hewlett  Packard  Model  44701A  Digital  Integrating  Voltmeter 

2.  Hewlett  Packard  Model  44705A  20  channel  relay  multiplexer  for  measuring  re¬ 
quired  voltage  drops. 

3.  Hewlett  Packard  Model  447 13A  24  channel  relay  multiplexer  with  electronic  cold 
junction  for  direct  conversion  of  thermal  couple  emf  to  temperature. 

The  HP300  desktop  computer  was  used  to  control  the  HP3852A  Data  acquisition  sys¬ 
tem.  Additionally,  through  specifically  written  programs,  the  computer  can  accomplish 
the  required  data  reductions,  (See  Appendix  D). 
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Figure  4.  Side  View  of  Experimental  Apparatus 


II 


Figure  5.  Rear  View  of  Experimental 
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III.  EXPERIMENTAL  PROCEDURE 


A.  PREPARATIONS 

Careful  preparations  were  required  to  ensure  that  the  data  obtained  was  of  high 
quality.  The  most  notable  requirement  being  the  necessity  of  thoughly  degassing  the 
FC-72  liquid,  [Ref,  p.  3]. 

1.  Normal  Warm-Up  Procedure 

In  carrying  out  the  experimental  runs  the  following  procedures  were  followed. 

1.  Power-up  the  HP3852A  Data  Acquisition  Unit  to  start  the  warm-up  period  of  one 
hour  minimum.  [Ref.  13:  p.  586] 

2.  Ensure  the  liquid  level  is  above  Wire  5,  but  below  the  condensate  drip  tray. 

3.  Set  power  supply  of  thermoelectric  coolers  to  1.3  amps  and  2.0  volts. 

4.  Ensure  the  experimental  chamber's  insulation  jacket  is  in  place  and  set  power  sup¬ 
ply  to  strip  heaters  to  approximately  1.8  amps  and  20  volts. 

5.  Verify  the  0.25  inch  diameter  vent  line  is  unobstructed. 

6.  Power  on  the  HP300  series  desktop  computer,  and  load  data  acquisition  program 
as  contained  in  Appendix  D. 

7.  Once  fluid  has  reached  saturation  and  is  boiling  vigourisly,  hold  for  a  minimum  of 
one  hour,  typically  two  hours  was  achieved,  to  permit  adequate  degassing  of  fluid. 

In  actual  practice,  the  time  required  to  raise  the  bulk  temperature  of  the  fluid  to  satu¬ 
ration  temperature  of  56  °C  was  nominally  two  hours. 

2.  Post  Warm-up  Procedure 

Once  the  required  degassing  has  been  accomplished,  power  is  temporarily  se¬ 
cured  to  the  strip  heaters.  Once  power  is  secured  the  forward  strip  heater  is  electrically 
isolated.  The  rear  strip  heater  is  powered  at  0.9  amps  and  20  volts,  and  the  thermo 
electric  coolers  are  secured. 

The  power  setting  of  the  strip  heaters  permits  vigorous  boiling  from  it's  surface. 
This  serves  to  keep  the  fluid  degassed  and  the  bulk  fluid  temperature  near  saturation  and 
near  isothermal.  By  keeping  the  fluid  level  below  the  insert  board's  upper  edge,  the  fluid 
forward  of  the  insert  board  remains  undisturbed  by  the  boiling  action  occurring  in  the 
rear  chamber. 
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B.  DATA  ACQUISITION  PROCEDURE 

Upon  completion  of  the  previous  preparatory  procedures,  the  experimental  runs  are 
ready  to  begin.  The  procedure  followed  for  an  experimental  run  is  as  follows: 

1.  Plug  power  supply  into  power  receptacles  for  wire  1.  Power  on  and  ensure  current 
limiter  is  set  to  max  value  of  1.0  amps. 

2.  Plug  HP6214C  power  supply  into  power  receptacles  of  desired  constant  heat  flux 
wire. 

3.  Set  power  supply  of  step  2  to  proper  output  to  give  desired  heat  flux  from  the  wire. 

4.  Start  running  the  Data  Acquisition  Program,  which  is  programmed  to  measure  the 
following  data  from  the  desired  wires  at  two  minute  intervals. 

•  Precision  Resistor  Voltage  drop,  up  to  five  channels. 

•  Platinum  Wire  Voltage  drop,  up  to  five  channels. 

•  Four  bulk  fluid  temperature  measurements 

•  Two  condenser  surface  temperature  measurements. 

•  One  ambient  air  temperature  measurement. 

5.  After  each  data  interval,  increase  the  voltage  to  the  varying  heat  flux  wire  by  0.2 
Volt  increments.  This  was  continued  until  the  on  line  power  supply  was  at  2.3  volts 
of  outputs,  the  increments  were  then  reduced  to  0.1  volts  until  nucleation  occurred. 
The  increments  were  then  returned  to  0.2  volts  until  a  heat  flux  of  approximately 
200,000  watts/m^  were  achieved.  At  any  higher  heat  flux  the  wire  will  go  into  the 
film  boiling  regime. 

6.  Reduce  the  power  to  the  varying  heat  flux  wire  at  0.3  Volts  increment  to  bring  the 
wire  back  to  a  no  power  condition. 

7.  Repeat  steps  3  thru  6  for  all  other  wire  and  heat  flux  configurations. 


C.  DATA  REDUCTION 


Data  that  was  obtained  from  the  actual  experimental  apparatus  were  two  voltage 
measurements  per  active  platinum  wire  channel,  and  six  temperature  measurements. 
The  process  used  for  data  reduction  is  a  combination  of  computer  programs  and 
graphical  outputs  for  numerical  presentation  and  visual  presentations  respectively. 

1.  Numerical  Data 

The  numerical  data  values  required  to  be  processed  are: 


1.  Fluid  Bulk  Temperature 


(T.  +  Tj  +  Ta  +  TJ 
^bulk-  4^0  t  '-1 


(3.1) 


and  where,  T,  thru  T4  are  the  thermocouple  measurements  of  the  bulk  fluid. 
2.  Current  flow  through  the  Platinum  Wire,  where 
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(Amps) 


(3.2) 


^wire ' 


'Res 


R 


2a 


3.  Platinum  Wire's  electrical  resistance  value,  where 

Rp/iwrc  =  (Ohms) 

^wlre 

4.  Platinum  Wire  surface  temperature,  where 

Tw  =  ^^Ptmre  +  ToC’C) 

where, 

T„,=  the  actual  wires  surface  temperature 
e  is  the  slope  of  the  calibration  curve 
To  is  the  calibration  curves  y-ax*s  intercept 

5.  Fluid  Film  Temperature,  where 

T  _ 

^Jtlm  2.0 


6.  Dissipated  Heat  Flux  of  platinum  wire,  where 


Iw/re  ^  ^PtWlre  ,  Watts  ^ 
1  •)  > 


Area 


■Surf 


m 


where,  Area5„,/=  n  x  lengthy, „  x 


(3.3) 

(3.4) 


(3.5) 


(3.6) 


The  above  numerical  data  were  all  derived  from  the  raw  data,  and  previously 
determined  calibration  curves.  The  following  table  provides  physical  dimensions  and 
temperature  calibration  coefficients  of  the  platinum  wires. 


Table  2.  PLATINUM  WIRE  PHYSICAL  PARAMETERS  FOR  EQUATION 


(3.4)  AND  EQUATION  (3.5) 


Wire 

e 

t; 

length 

diameter 

'’C/Q 

"C 

mm 

mm 

1 

53.902 

-254.470 

99.22 

0.05 

2 

53.666 

-253.009 

99.22 

0.05 

3 

53.141 

-252.729 

100.01 

0.05 

4 

53.484 

-252.848 

98.43 

0.05 

5 

53.086 

-260.049 

100.01 

0.05 

To  see  in  detail  the  methodology  used  to  obtain  the  numerical  values  presented 
in  the  preceding  table,  proceed  to  Appendix  A  'Platinum  Wire  Calibration'. 
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2.  Fluid  Property  Data 

The  Fluid  Properties  to  be  determined  are  as  follows: 

1.  Thermal  Conductivity  of  fluid,  as  obtained  from  Figure  5  of  the  3M  Corporation 
Flourinert  Product  Manual  [Ref.  3:  p.  15]  have  been  determined  to  be, 


0.6033 -0.001 15  xTyi„„ 

io! 


( 


Watts  .. 
°C  ^ 


(3.7) 


2.  Tliermal  Expansion  Coefficient,  from  table  4B  in  the  Product  Manual  [Ref.  3:  p.  10] 
is  given  as 


0.00261  ,  1  . 

^  1.740-0.00251  ^  “C  ^ 


(3.8) 


3.  Kinematic  Viscosity  of  the  Fluid  ,  as  obtained  from  Figure  (3)  in  the  Product  Man¬ 
ual  [Ref  3:.  p.l3|  and  determining  an  exponential  curve  fit  yields 


v=  1.203952F-8  x  exp(l, 058.4109° 


K 

T,„,,  +  273.15°C 


(3.9) 


4.  Specific  Heat  of  the  Fluid  for  all  the  Flourinert  Fluids  are  obtained  from  Figure  (4) 
of  the  Product  Manual  [Ref  3  :  p.  14]  is, 

Cp  =  (0.24111  -F  ^J0mE-4x  x  4186(-^^)  (3.10) 


5.  Liquid  Density  of  FC-72,  from  expression  in  table  4B  and  constants  from  table  4C 
in  the  Product  Manual  [Ref  3:  p  10]  is  given  by 


p  =  (1.740  -  0.00261  X  x  1000(  ) 


m 


6.  Thermal  Diffusivity  of  the  liquid  is  given  by. 


(3.11) 


a  = 


pxCp 


(3.12) 


7.  Prandlt  Number  of  the  Fluid  is  given  by 


8.  Ralyeigh  Number  of  Fluid  based  on  wire  diameter  is  given  by, 


Ra  =  - 


gx/?x(T^,-Tt„,*)xD 


V  X  a 


(3.13) 


(3.14) 


9.  Correlation  for  the  Nusselt  Number  was  obtained  by  using  the  correlation  of  Kuehn 
and  Goldstein  [Ref  10;  p.  1128]  given  by 


Nu 


D- 


lnCl-1- 


[(0.518  X  Kaf  X  [1  +  (-Mp-)3/5]-s/i2)i5  +  (q.i  x 


] 
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10.  Theoretical  Heat  Flux  as  obtained  from  using  the  value  obtained  from  equation 
(3.15). 


k  X  Nuo  X  ,  Watts  ^ 

q - 5 - (— r-) 


m 


(3.16) 


With  all  the  required  numerical  data  and  fluid  properties  obtained,  it  is  now 
necessary  to  present  the  numerical  information  in  a  graphical  representation  to  see  the 
relevancy  of  the  results  obtained.  This  will  become  apparent  with  the  presentation  of 
results  in  Chapter  IV  'Results  and  Discussion'. 
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IV.  RESULTS  AND  DISCUSSION 


The  data  obtained  can  be  grouped  for  discussion  as  follows: 

1.  Individually  powered  wires 

Saturated  Conditions 
Sub-Cooled  Conditions 

2.  Bubble  Pumping  Effect  Runs 

Saturated  Conditions 
Sub-Cooled  Conditions 

A.  INDIVIDUALLY  POWERED  WIRES 

Before  analyzing  the  pumping  effect  of  the  bubbles  on  the  upper  wires,  the  individ¬ 
ual  wire  boiling  experiments  will  be  investigated.  All  of  the  individual  powered  wire  runs 
shown  here  were  made  at  or  near  saturated  fluid  conditions  with  two  exceptions.  The 
following  table  will  provide  a  listing  of  single  wire  runs. 


Table  3.  INDIVIDUAL  WIRE  RUNS 


Figure 

Wire 

Figure  10 

1 

55,5"  C 

Figure  11 

1  -  several  runs  overlayed 

At  or  near  saturation 

Figure  12 

1 

38"  C 

Figure  13 

1 

23.0"  C 

Figure  14 

2 

56.0"  C 

Figure  15 

3 

56.0"  C 

Figure  16 

4 

56.0"  C 

Figure  17 

all  the  above 

see  individual  wires 

In  all  cases  conditions  were  under  atmospheric  pressure  with  0.05mm  diameter  platinum 
wires. 

1.  Natural  Convection 

Figure  10  is  a  typical  boiling  curve  for  FC-72  with  the  Kuehn  and  Goldstein 
natural  convection  region  shown  along  with  the  experimental  points.  The  nucleate 
boiling  zone  starts  at  boiling  incipience,  which  is  evident  from  the  sudden  drop  of  wall 
temperature,  a  characteristic  of  a  highly  wetting  fluid.  Figure  1 1'is  a  compilation  of  wire 
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1  data,  at  or  near  saturation,  obtained  over  a  period  of  time.  As  can  be  seen,  with  the 
exception  of  the  boiling  incipience  point,  the  data  is  virtually  super  imposable  over  each 
other.  This  reproducability  of  data  over  time  indicates  no  significant  aging  has  occured 
on  wire  1,  and  more  importantly  no  change  in  calibration. 

In  fact  in  the  natural  convection  region  Figure  10  thru  Figure  16  arc  ah  essen¬ 
tially  superimposable  over  each  other.  Figure  17  is  an  overlay  of  all  the  previous  fig¬ 
ures.  Therefore,  it  appears  that  whether  the  fluid  is  Saturated  or  Sub-Cooled,  the  fluid 
bulk  temperature  appears  to  have  little  effect  on  the  natural  convection  heat  transfer 
rate. 

There  was  generally  very  good  agreement  between  the  natural  convection  data 
and  the  Kuehn  and  Goldstein  correlation  [Ref  10]  (equation  3.15).  The  only  deviation 
existing  between  the  figures  being  the  total  temperature  excursion  to  incipience,  which 
ranged  from  2T  C  to  38°  C. 

2.  Nucleate  Boiling 

As  previously  stated,  the  nucleate  boiling  zone  starts  at  boiling  incipience.  Us¬ 
ing  wire  1  as  a  standard  and  Figure  lOas  base  line  boiling  characteristic  data  a  compar¬ 
ison  to  other  vires  at  saturation  condition,  and  a  comparison  to  subcooling  effects  on 
wire's  1  boiling  can  be  made.  As  is  evident  from  Figure  17,  the  boiling  characteristic 
of  all  the  wires  at  saturation  condition  tend  to  overlap  each  other.  The  differences  aris¬ 
ing  tending  to  be  caused  by  individual  wire  placement  in  the  fluid  .  However,  the  effect 
of  subcooling  does  result  in  a  decrease  in  the  boiling  curve  overshoot  at  point  of 
incipience.  This  decrease  in  the  boiling  curve  overshoot  results  in  a  net  decrease  in  the 
superheat  required  to  initiate  nucleate  boiling. 

3.  Departure  from  Nucleate  Boiling 

The  Critical  Heat  Flux  (CHF)  for  FC-72  at  saturated  condition  was  calculated 
by  Danielson  et  all  [Ref  9]  to  be  approximately  203,000  Watts/m^  Experimental  value 
for  CHF  ranged  from  198,000  Watts/m^  to  220,000  Watts/mS  a  spread  of  -2.5%  to 
8.4%.  A  CHF  for  a  subcooled  condition  of  Ti„tt=38°  C  was  found  experimentally  to 
be  approximately  400,000  watts/m^  A  CHF  for  T*„,*=  23  °  C  was  not  obtained  during 
the  experimental  runs. 
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HEAT  FLUX  (WATTS/M2) 


WIRE  1  BOlUNG  CURVE 


NEAR  SATURATION  =  S5.5°Celisus 


Figure  10.  Wirel  boiling  cune  Ti„tt  =  55.5°C 
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WIRE  1  BOlUNG  CURVES 


Figure  11.  Wire  1  Boiling  Curves  over  time 
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HEAT  FLUX  (WATTS/M2) 
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HEAT  FLUX  (WATTS/M2) 


WIRE  I  eOIUMG  CURVE 


Figure  13.  Wire  1  Boiling  Curve  Tj.,*  =  23°C 
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HEAT  FLUX  (WATTS/M2) 


WIRE  3  BOILING  CURVE 


S  ATURAT  =  56.0”  Celbus 


Figure  15.  Wire  3  Boiling  Curve  Tj.,*  =  56‘’C 
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HEAT  FLUX  (WATTS/M2) 


WIRE  4  BOlUNG  CURVE 


SATURATED  =  56.0°  Celisus 


Figure  16.  Wire  4  boiling  curve  T^^a  =  56'’C 
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HEAT  FLUX  (WAnS/M2) 


BOILING  CURVES  COMPILED 
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Figure  17.  Overlay  of  previous  curves 
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B.  BUBBLE  PUMPING  EFFECT  ON  HEAT  TRANSFER  RATE 


In  the  upcoming  discussion  a  brief  explanation  of  the  figures  are  required.  In  each 
case  for  these  experiments,  heater  wire  number  I  was  used  to  generate  the  bubbles  for 
examining  the  pumping  effects  upon  the  upper  wires.  The  figure  insert  is  used  to  show 
which  wires  are  being  utilized.  Only  select  data  points  were  plotted  for  the  passive 
(constant  heat  flux  wire)  wire,  and  those  points  correspond  in  time  to  the  identically 
numbered  data  point  on  wire  I's  boiling  curve.  The  data  points  were  chosen  since  they 
define  regions  of  wire  I's  boiling  curve  as  shown  in  the  following  table. 


Table  4.  DATA  POINT  DESCRIPTION 


5  point  figure 

6  point  figure 

Region  Description 

1 

1 

Natural  Convection  Start 

2 

2 

Point  of  Incipience 

3 

Partial  Boiling  on  wire 

3 

4 

Full  nucleate  Boiling 

4 

5 

Max  Nucleate  Boiling 

5 

6 

End  of  run 

1.  Passive  Wire  at  Low  Heat  Flux 

In  these  experimental  runs  the  heat  flux  on  the  passive  wire  was  set  to  approx¬ 
imately  8500  watts/m^  At  this  low  heat  flux  it  was  expected  that  little  to  no  active 
nucleation  will  occur  on  the  passive  wire,  even  after  wire  1  is  in  nucleate  boiling. 
Therefore,  all  the  enhancement  of  the  upper  wire  will  be  due  to  an  addition  of  a  force 
convection  term  to  the  heat  transfer  rate. 

In  comparing  Figure  18  to  Figure  19the  following  can  be  inferred: 

1.  Total  reduction  of  wall  temperature  on  the  passive  wires  are  comparable.  Wire  2 
having  a  3.7°C  and  Wire  4  having  a  3.0®  C  drop  in  wall  temperature. 

2.  A  definite  effective  distance  exist  over  which  wire's  1  influence  is  imparted  upon  the 
upper  wires,  for  example: 

Wire  2  and  wire  4  showed  a  very  small  change  in  wall  temperature  at  point  2  (point 
of  incipience). 

After  incipient  boiling  has  occurred  on  wire  1,  wire  2  showed  a  greater  drop  in  wall 
surface  temperature  then  wire  4.  Approx.  2.5°  C  drop  for  wire  2  and  0.1°  C  drop 
for  wire  4. 
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3.  At  point  4,  the  point  of  maximum  temperature  in  nucleate  boiling,  both  wire  2  and 
wire  4  showed  a  comparable  surface  temperature,  approximately  T  C  above  satu¬ 
ration. 

As  can  be  seen  from  the  figures,  some  nucleation  sites  have  been  activated  on  wire  2  and 
wire  4,  with  point  5  showing  a  lower  wall  temperature  then  point  1. 

2.  Passive  Wire  at  Medium  Heat  Flux 

In  these  experimental  runs  the  heat  flux  on  the  passive  wire  was  set  to  approx¬ 
imately  20,500  watts/ml  At  this  heat  flux  it  was  expected  that  active  nucleation  will 
occur  on  the  passive  wire  after  wire  1  is  nucleate  in  boiling. 

In  comparing  Figure  20  to  Figure  21 
the  following  can  be  inferred: 

1.  The  total  reduction  of  wall  temperature  on  the  passive  wire  ,  points  1  to  points  4, 
are  comparable  at  7.1®  C  for  wire  2,  and  6.3®  C  for  wire  4. 

2.  The  spatial  effect  is  evident  on  the  influence  of  wire  1  on  the  upper  wires,  for  ex¬ 
ample: 

Between  point  1  and  point  2,  wire  2  shows  a  greater  reduction  in  wall  surface  tem¬ 
perature  than  is  evident  in  wire  4. 

At  point  3,  wire  2  shows  a  much  greater  drop  in  surface  temperature  than  wire  4: 
4®  C  compared  to  2®  C.  The  difference  being  the  amount  of 

3.  At  point  4,  maximum  nucleate  boiling  from  wire  1,  both  wire  2  and  wire  4  show 
comparable  values  for  surface  temperature,  approximately 

1.3®  C. 

As  can  be  seen  from  the  figures,  wire  1  has  activated  previously  inactive  nucleation  sites 
on  wire  2  and  wire  4.  This  is  evident  from  point  5  being  at  a  lower  surface  temperature 
than  point  3. 

3.  Passive  Wire  at  High  Heat  Flux 

In  these  experimental  runs  the  heat  flux  on  the  passive  wire  was  set  to  40,000 
Watts/m^  on  wire2  and  35,000  Watts/m^  on  Wire  4.  At  this  heat  flux  it  was  expected 
that  a  large  number  of  inactive  nucleation  sites  will  be  activated  after  wire  1  is  nucleate 
in  boiling. 

In  comparing  Figure  22  and  Figure  23  it  is  evident  that  the  information  ob¬ 
tainable  from  the  figures  is  comparable  to  the  information  for  the  medium  heat  flux  case. 
However,  the  wall  temperature  reductions  on  the  passive  wires  are  much  greater  as 
compared  to  the  Medium  Heat  Flux  case,  in  particular  the  following: 

I.  In  comparing  the  surface  temperature  reduction  from  point  1  to  point  2,  wire  2  has 
dropped  2.5®  C  while  wire  4  had  almost  a  zero  drop  of  temperature.  This  appears 
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to  be  due  to  the  relative  locations  of  the  passive  wires  with  respect  to  wire  1,  and 
is  especially  apparent  at  the  higher  heat  fluxes. 

2.  Many  inactive  nucleation  sites  have  been  activated  by  wire  I's  nucleation  bubbles. 
As  can  be  seen  by  comparing  point  1  to  point  5(6),  a  difference  on  Wire  2  of  13® 
C  and  on  Wire  4  of  7°  C  is  apparent. 

4.  Wire  I  Passive  with  Wire  2  varying 

In  these  experimental  runs,  wire  1  was  set  to  150,000  Watts/m^  so  it  was  boiling 
vigorously,  wire  2  was  actively  varied  over  a  large  heat  flux  range  as  can  be  seen  from 
Figure  24.  On  Figure  24  a  previously  obtained  boiling  curve  of  Wire  2  was  overlayed 
onto  it.  Upon  comparing  the  two  curves  the  most  striking  feature  is  the  disappearance 
of  the  Hysterises  Loop  as  is  evident  in  the  overlayed  boihng  curve  data.  In  this  ar¬ 
rangement,  the  boiling  curve  overshoot  has  been  removed  as  a  design  problem,  bubbles 
from  Wire  1  has  removed  a  major  design  consideration. 

5.  Subcooled 

In  these  experiments,  the  bulk  fluid  was  degassed  and  then  subcooled  to  a  tem¬ 
perature  of  38°  C.  In  these  experimental  runs,  wire  2  was  set  at  35,000  watts/m^  and 
53,000  watts/mS  as  shown  in  Figure  25  and  Figure  26  respectively.  In  subcooling,  no 
significant  effects  were  notable  beyond  those  seen  in  the  saturated  runs.  The  only  effect 
that  was  apparent  was  the  ability  to  run  wire  2  at  a  greater  heat  flux  with  out  incipent 
boiling  occuring  on  the  wire. 
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HEAT  FLUX  (WAnS/M2) 
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HEAT  FLUX  (WATTS/M2) 
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HEAT  FLUX  (WAnS/M2) 
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HEAT  FLUX  (WATTS/M2) 


WIRE  4  AT  20,500  WArT/M2 


NEAR  SATURATION  T,,,,  =  SS.rCelisus 


Figure  21.  Wire  4  at  20,500  W/  T*^,*  =  55.6X 
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HEAT  FLUX  (WATTS/M2) 


HEAT  FLUX  (WAnS/M2) 


24.  Wire  1  passive,  Wire  2  active  Tj.,*  =  56°C 


HEAT  FLUX  (WATTS/M2) 
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HEAT  FLUX  (WATTS/M2) 


WIRE  2  AT  55.C00  WATT/W2 


SUBCOOLED  =  3S° Celisus 


Figure  26.  Wire  2  at  53,000  W/  =  SS'C 
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V.  CONCLUSIONS 


The  study  conducted  in  this  thesis  was  to  investigate  the  efTects  of  the  wake  from  a 
bubble  generating  wire  on  wires  located  dc'”nstream  within  the  wake.  Experimental 
runs  were  conducted  at  saturated  condition  and  subcooled  conditions  in  which  the  fol¬ 
lowing  were  found: 

A.  SINGLE  WIRE  BOILING 

1.  The  effects  of  subcooling  on  the  natural  convective  heat  transfer  from  the  wire  was 
negligible. 

2.  The  effects  of  subcooling  tended  to  increase  the  heat  flux  required  for  incipient 
boiling  to  initiate. 

3.  With  subcooling,  boiling  hystersis  is  still  observed,  but  the  hystersis  magnitude  de¬ 
creases  as  subcooling  increases. 

4.  Subcooling  results  in  a  reduction  in  wall  temperature  after  the  start  of  nucleate 
boiling  as  compared  to  saturated  cases. 

B.  PASSIVE  WIRE  ANALYSIS 

Heat  Transfer  enhancement  effects  from  Wire  i  upon  the  downstream  wires  has 
many  factors  within  it.  The  following  factors  have  been  identified: 

1.  An  enhancement  to  convection  resulting  from  buoyancy  induced  fluid  flow  off  of 
wire  1  was  seen  to  exist.  This  effect  is  very  spatially  dependent  with  the  effects 
most  pronounced  on  wire  2  and  negligible  on  wire  4.  This  corresponds  to  point  1 
and  point  2  on  the  figures. 

2.  The  sudden  drop  of  wall  temperature  between  point  2  and  point  3/4  was  seen  to 
be  the  results  of  the  following: 

Nucleation  bubbles  from  wire  1  activating  the  inactive  nucleation  sites  on  the  upper 
wires  resulting  in  boiling  heat  transfer  off  of  the  upper  wires. 

An  additional  force  convection  enhancement  due  to  the  nucleation  bubbles  from 
wire  1  breaking  up  the  thermal  boundary  layer  about  the  upper  wires  is  believed  to 
exist. 

3.  The  drop  in  temperature  on  the  upper  wires  corresponding  with  points  3/4  to  point 
4/5  appears  to  be  due  purely  to  an  increase  in  force  convection  due  to  the  high 
density  of  bubbles  off  of  wire  1. 

Additional  conclusions  reached  are. 

1.  A  constant  source  of  bubble  generation  from  below  is  an  effective  way  to  eliminate 
the  Boiling  Curve  Overshoot  and  Hysteresis  loop  associated  with  dielectric  fluids. 


44 


2.  At  low  heat  fluxes,  proximity  to  Wire  1  is  important  for  heat  transfer  enhancement 
to  occur.  This  dependency  reduces  as  heat  flux  increase  on  the  passive  wire. 
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APPENDIX  A.  CALIBRATION 


A.  PLATINUM  WIRE 

1.  Surface  Micro-Structure 

Prior  to  the  experiments,  a  length  of  0.05mm  platinum  wire  was  analyzed  using 
a  Scanning  Electron  Microscope  (SEM).  The  length  of  wire  analyzed  was  not  used  in 
any  experiments,  but  it  did  come  from  the  same  spool  of  wire  as  the  experimental  lengths 
of  wire,  therefore,  it  should  be  a  representative  sample.  In  analyzing  this  sample  the 
SEM  brought  out  several  important  features: 

1.  From  Figure  27  a  wire  diameter  of  0.052mm  can  be  determined.  Also,  the  surface 
is  relatively  unmarked  at  this  magnification  of  1490x;  little  scoring  on  the  wire  and 
a  low  density  of  surface  pits  are  evident. 

2.  In  Figure  28,  an  increase  in  magnification  to  4740x  within  the  field  of  Figure  27, 
a  more  detailed  view  of  the  larger  surface  defects  are  evident.  Maximum  defect  size 
is  approx  2/xm  in  length  and  0.5/xm  in  width.  A  surface  defect  of  that  size  would 
be  too  large  to  be  a  possible  nucleation  site  ,  according  to  the  calculations  done  in 
Ref.  4  ofO.l/xm. 

3.  In  Figure  29,  a  further  increase  in  magnification  to  13,200x  within  the  field  of 
Figure  27,  pits  of  approximately  O.S^m  max  down  to  0.3^m  are  seen.  It  will  be 
these  surface  defects  that  boiling  will  initiate  from.  The  density  of  these  pits  along 
the  length  of  the  wire  analyzed  was  not  great,  this  implies  that  the  incipient  super¬ 
heat  wUl  be  large. 


Figure  27.  Surface  micro  geometry:  Sample  Platinum  Wire  (1490x) 
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Figure  28.  Surface  micro  geometry:  Sample  Platinum  Wire  (4740x) 


Figure  29.  Surface  micro  geometry:  Sample  Platinum  Wire  (13,200x) 


2.  Desirability  for  Using  Platinum 

Platinum  was  chosen  for  several  reasons: 

1.  It  has  a  highly  linear  response  of  Temperature  vs  Resistance. 

2.  It  has  a  large  AR  for  a  given  AT. 

In  each  aspect,  platinum  performed  very  well.  As  can  be  seen  in  the  calibration  curves 
of  Figure  31  thru  Figure  40  response  is  linear,  with  a  0.1  £2  change  in  resistance  for  ev¬ 
ery  5°  C  change  in  temperature. 

B.  BACKGROUND 

An  important  step  in  minimizing  experimental  data  uncertainty  is  an  accurate  cali¬ 
bration  of  the  platinum  wires  used.  In  obtaining  the  calibration  curves  for  the  wires, 
three  calibration  runs  were  performed. 

1.  An  initial  calibration  in  a  heated  calibration  bath  with  a  platinum  resistance 
thermometer. 

2.  An  in  place  calibration  of  wires  in  the  test  chamber,  with  the  four  installed 
thermocouples  used  for  temperature  measurements,  was  conducted  to  ensure  the 
calibration  curves  from  the  calibration  bath  were  not  altered  from  handling  the  wire 
insert  board. 

3.  A  final  in  place  calibration  of  wires  in  the  test  chamber, 

with  the  four  installed  thermocouples  used  for  temperature  measurements,  was 
conducted  to  ensure  the  calibration  curves  used  for  data  acquisition  were  not  al¬ 
tered  due  to  any  reason. 

The  second  calibration  run  was  required  to  ensure  no  distortion  occurred  to  the  wires 
when  the  insert  board,  on  which  the  wires  are  mounted,  was  removed  from  the  calibrated 
bath  and  installed  into  the  experimental  test  chamber. 

C.  INITIAL  CALIBRATION 

The  initial  calibration  was  conducted  using  the  following  equipment. 

1.  Rosemount  Engineering  Co.  Model  913A  calibration  bath  with  Etheylene  Glycol 
as  the  working  fluid. 

2.  Rosemount  Engineering  Co.  Model  923B  power  supply,  controlling  calibration 
bath  temperature. 

3.  Rosemount  Engineering  Co.  Model  920A  commutating  bridge,  used  to  measure  an 
resistance  to  0.001£2  accuracy. 

4.  Rosemount  Engineering  Co.  Model  162C  S/N  985  Platinum  Temperature  Probe. 
Used  to  measure  calibration  bath  temperature  through  the  commutating  bridge. 

5.  Hewlett  Packard  Data  Acquisition  System.  Used  to  measure  platinum  wire  resist¬ 
ance  using  a  4-wire  ohm  measurement  technique. 
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1.  Calibration  Heating  Bath 

The  present  investigation  was  conducted  using  an  externally  heated  calibration 
bath  with  etheylene  glycol  as  the  working  fluid.  This  fluid  was  circulated  by  a  centrifugal 
pump  to  ensure  a  uniform  bulk  temperature  throughout  the  bath.  Item  (1)  and  item  (2) 
are  utilized  . 

2.  Temperature  Measurements 

Temperature  measurements  for  calibration  were  conducted  using  a  high  sensi¬ 
tivity  commutating  bridge  to  measure  the  resistance  of  a  platinum  resistance  temper¬ 
ature  probe.  This  platinum  temperature  probe  had  been  previously  calibrated  to  a 
standard  that  is  traceable  to  the  National  Bureau  of  Standards. 

3.  Resistance  Measurements 

In  calibrating  the  platinum  wire  the  resistance  of  the  wire  was  measured  with  the 
same  Data  Acquisition  Unit  Components  as  would  be  used  in  the  experimental  runs. 
This  was  to  ensure  the  uncertainty  of  measured  resistance  would  be  the  same  as  between 
the  calibration  run  and  experimental  runs. 

The  Data  Acquisition  Unit  Component  used  are; 

1.  HP3852A  Control  Unit 

2.  HP4470iA  Integrating  Voltmeter  card 

3.  HP44705A  20  channel  Multiplexer  card 

4.  HP300  series  desk  top  computer 

5.  HP44713A  24  channel  multiplexer  with  electronic  cold  junction  for  thermo  couple 
measurements. 

4.  Calibration  Procedure 

By  using  the  previously  described  components  for  measuring  temperature  and 
resistance,  the  following  procedure  was  used  to  obtain  the  system's  calibration  data: 

1.  Set  calibration  bath  controller  for  35®  C 

2.  Allow  bath  temperature  to  stabilize. 

3.  Carefully  place  Insert  Board  into  bath  liquid.  EXTREME  care  must  be  taken  to 
ensure  that  none  of  the  wires  are  deformed  or  destroyed. 

4.  Allow  insert  board  to  come  to  thermal  equilibrium  with  the  bath, 
this  takes  approximately  five  minutes. 

5.  Take  temperature  reading  with  platinum  wire  temperature  probe  and  commutator 
bridge,  and  resistance  measurement  of  wires  with  data  acquisition  unit. 

6.  Remove  insert  board,  and  set  bath  temperature  controller  for  a  5®  C  increase. 

7.  Repeat  step  2  thru  6  until  an  90®  C  reading  has  been  achieved. 
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8.  Activate  External  Cooling  system. 

9.  With  temperature  control  setting  decreasing  at  5°  C  intervals,  repeat  step  2  thru  6 
until  a  35®  C  reading  has  been  achieved. 

D.  IN  PLACE  CALIBRATION 

To  ensure  that  the  calibration  curve  of  platinum  wire  1  thru  5  are  reliable,  an  in 
place  calibration  run  was  performed  on  the  wires.  This  calibration  procedure  was  per¬ 
formed  before  any  experimental  runs  were  conducted,  and  after  the  last  experimental  run 
was  completed  to  detect  any  changes  in  calibration  due  to  wire  deformation  or  aging. 

The  equipment  utilized  for  the  inplace  calibration  was  the  same  Data  Acquisition 
equipment  used  for  the  resistance  measurements  of  the  Bath  CaUbration  runs  with  two 
additions. 

1.  HP44713A  24  Channel  Multiplexer  with  electronic  cold  junction. 

2.  Four  type  T  thermocouples. 

The  thermocouples  are  now  utilized,  instead  of  the  platinum  resistance  temperature 
probe,  for  measuring  fluid  bulk  temperature.  The  uncertainity  associated  with  the 
thermocouples  is  0.4®C  as  compared  to  the  platinum  resistance  temperature  probe  used 
as  a  standard  Figure  30  shows  the  thermocouples  uncertainity  compared  to  the 
platinum  resistance  temperature  probe. 

I.  Inplace  Calibration  Procedure. 

The  following  procedure  is  used  for  inplace  calibration  of  the  platinum  wires: 

1.  Ensure  Data  Acquisition  system  is  wired  for  4-wire  ohm  measurements. 

2.  Turn  on  Data  Acquisition  system  and  allow  one  hour  of  warm-up  time. 

3.  Load  Calibration  Acquisition  Program  into  the  HP-300  series  Computer.  See  Ap¬ 
pendix  D  for  program  listing. 

4.  Set  power  supply  to  bulk  heaters  to  20volts  and  1.8  amps. 

5.  Start  Acquisition  program.  The  acquisition  program  will  measure  the  resistance 
of  each  platinum  wire,  and  measure  the  temperature  output  of  each  thermocouple 
at  a  10  minute  interval. 

6.  Allow  acquisition  program  to  run  until  bulk  fluid  has  reached  Saturation  temper¬ 
ature  of  56°  C. 

E.  CALIBRATION  RESULTS 
1.  Platinum  Wire  1 

From  Figure  31  the  importance  of  conducting  an  in  place  calibration  prior  to 
any  experimental  runs  becomes  obvious,  wire  I's  Bath  calibration  data  has  been  con- 
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siderably  altered  wliile  it  was  being  inserted  into  the  test  chamber.  The  initial  in  place 
calibration  data  was  used  for  the  experimental  runs  calibration  curve.  The  resulting 
calibration  was: 


T,  =  53.902  X  R,  -  254.470  (‘’C) 

From  Figure  32  it  can  be  seen  that  the  post  calibration  curve  matches  well  with 
the  data  acquisition  curve. 

2.  Platinum  Wire  2 

From  Figure  33,  it  is  apparent  that  the  initial  in  place  calibration  curve  fairs  in 
well  with  the  calibration  bath's  decreasing  data.  The  initial  in  place  calibration  data  was 
used  for  the  experimental  runs  calibration  curve.  The  resulting  calibration  was: 

Tj  =  53.666  X  Rj  -  253.009  (°C) 

From  Figure  34  it  can  be  seen  that  the  Post  run  calibration  matches  well  with 
the  Data  Acquisition  Calibration  curve. 

3.  Platinum  Wire  3 

From  Figure  35  it  is  apparent  that  wire  3's  Bath  calibration  data  has  been  al¬ 
tered  between  the  calibration  bath  and  insertion  into  the  test  chamber.  The  initial  in 
place  calibration  data  was  used  for  the  experimental  runs  calibration  curve.  The  result¬ 
ing  calibration  was: 


T3  =  53.141  X  R3- 252.729  (X) 

From  Figure  36  it  can  be  seen  that  the  Post  run  calibration  matches  well  with 
the  Data  Acquisition  Calibration  curve. 

4.  Platinum  Wire  4 

From  Figure  37,  it  is  apparent  that  the  initial  in  place  calibration  curve  fairs  in 
well  with  the  calibration  bath's  decreasing  data.  The  initial  in  place  calibration  data  was 
used  for  the  experimental  runs  calibration  curve.  The  resulting  calibration  was: 

T4  =  53.484  X  R4  -  252.848  (X) 

Form  Figure  38  it  can  be  seen  that  the  post  run  calibration  matches  well  with 
the  data  acquisition  calibration  curve. 
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5.  Platinum  Wire  5 

From  Figure  39  it  is  clear  to  see  the  shift  in  calibration  curves  between  the  cal¬ 
ibration  bath  data  and  initial  inplace  calibration  data.  The  calibration  data  has  been 
altered  while  it  was  being  inserted  into  the  test  chamber  obviously.  The  initial  data 
calibration  was  used  for  the  experimental  run  calibration  curve.  The  resulting  cali¬ 
bration  was: 

Tj  =  53.085  X  R5  -  260.049  (X) 

From  Figure  40  it  can  be  further  seen  that  a  considerable  error  has  been  intro¬ 
duced  into  wire  5's  calibration,  an  error  of  3.7°  C.  This  places  an  uncertainity  value  onto 
wire's  5  temperature  daca  of  at  least  3.7°  C,  an  unacceptable  value. 
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THERMOCOUPLE  CALIBRATION 


Figure  30.  Thermocouple  Uncertainity:  Uncertainity  of  thermocouple's 
tempertaure  measurements  as  compared  to  a  platinum  resistance  tem¬ 
perature  probe. 
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PLATINUM  TEMPERATURE  STANDARD 


(swHo)  aoNvisisaa 


Figure  31.  Calibration  of  Platinum  Wire  1:  Calibration  curve  of  wire  I  of  Cali¬ 
bration  Bath  Temperature  Increasing,  Calibration  Bath  Temperature 
Decreasmg,  and  of  the  Initial  in  place  calibration  data. 
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WIRE  1  CALIBRATION 


1  - 

Figure  32.  Calibration  of  Platinum  Wire  1:  Calibration  curve  of  wire  I  of  Initial 


in  place  calibration  curve,  and  of  Final  in  place  calibration  curve 
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'I’EMPERATURE  (DEG  C) 


WIRE  2  CALIBRATION 
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Figure  33.  Calibration  of  Platinum  Wire  2:  Calibration  curve  of  wire  2  of  Cali¬ 
bration  Bath  Temperature  Increasing,  Calibration  Bath  Temperature 
Decreasing,  and  of  the  Initial  in  place  calibration  data. 
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TEMPERATURE  (DEG  C) 


WIRE  2  CALIBRATION 


Figure  34.  Calibration  of  Platinum  Wire  2:  Calibration  curve  of  wire  2  of  Initial 
in  place  calibration  curve,  and  of  Final  in  place  calibration  curve 
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TEMPERATURE  (DEG  C) 


WIRE  3  CALIBRATION 


Figure  35.  Calibration  of  Platinum  Wire  3:  Calibration  curve  of  wire  3  of  Cali¬ 
bration  Bath  Temperature  Increasing,  Calibration  Bath  Temperature 
Decreasing,  and  of  the  Initial  in  place  calibration  data. 
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TEMPERATURE  (DEG  C) 


WIRE  3  CALIBRATION 


in  place  calibration  curve,  and  of  Final  in  place  calibration  curve 
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WIRE  4  CALIBRATION 
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Figure  37.  Calibration  of  Platinum  Wire  4;  Calibration  curve  of  wire  4  of  Cali¬ 
bration  Bath  Temperature  Increasing,  Calibration  Bath  Temperature 
Decreasing,  and  of  the  Initial  in  place  calibration  data. 
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WIRE  4  CALIBRATION 


Figure  38.  Calibration  of  Platinum  Wire  4:  Calibration  curve  of  wire  4  of  Initial 


in  place  calibration  curve,  and  of  Final  in  place  calibration  curve 
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WIRE  5  CALIBRATION 


f  oil  UV  £4  •a  V  bli 

(swHo)  aoMvxsisaa 


Figure  39.  Calibration  of  Platinum  Wire  5:  Calibration  curve  of  wire  5  of  Cali¬ 
bration  Bath  Temperature  Increasing,  Calibration  Bath  Temperature 
Decreasing,  and  of  the  Initial  in  place  calibration  data. 
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WIRE  5  CALIBRATION 
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Figure  40.  Calibration  of  Platinum  Wire  5:  Calibration  curve  of  wire  5  of  Initial 
in  place  calibration  curve,  and  of  Final  in  place  calibration  curve 

F.  SUMMARY 

The  following  points  were  brought  out  during  the  calibration  process. 

1.  The  process  of  moving  the  platinum  wires  from  the  calibration  bath  and  installing 
them  into  the  test  chamber  resulted  in  a  shift  in  the  calibration  curves. 
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TEMPERATURE  (DEG  C) 


2.  The  most  probable  cause  of  all  the  calibration  shifts  were  the  lead  wires  in  the  back 
of  the  insert  board  being  pulled  and  adding  an  additional  tension  onto  the  platinum 
wires,  thereby  stretching  them.  The  fact  that  the  calibration  curves  cf  ’stently 
shifted  to  an  increase  in  resistance  supports  this  theory. 

3.  The  maximum  difference  observed  between  the  Platinum  Wire  Temperature  Probe 
and  the  Type  T  thermocouples  were  0.4°  C.  Therefore,  a  temperature  uncertainity 
for  the  in  place  calibration  will  be  set  at  0.4°  C. 

4.  The  differences  between  the  wire's  data  acquisition  calibration  curve  and  the  wire's 
final  calibration  curve  did  not  exceed  the  thermocouples  uncertainity  in  temper¬ 
ature. 
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APPENDIX  B.  SAMPLE  CALCULATIONS 


In  these  sample  calculations  the  output  value  of  wire  1  will  be  used. 

A.  DETERMINATION  OF  DERIVED  VALUES 

1.  Determination  of  the  Average  Bulk  Temperature 

The  average  bulk  temperature  of  the  liquid  was  obtained  by  the  arithmetic  av¬ 
erage  of  the  four  immersed  themocouple  outputs. 

(Ti  +  -H3  +  TJ 

^  4.0 

(55.835  +  55.884  +  55.988  +  55.932) 

4.0 

=  55.9098‘’C 

2.  Determination  of  Wire  Current 

The  current  flow  through  the  platinum  wire  is  obtained  by  measuring  the  volt¬ 
age  drop  across  a  precision  2  ohm  resistor. 

I  _ 

^PtWire  ~  n 

(0.74605VoIts) 

(l.99940hm) 

=  0.23810hms 

B.  DETERMINATION  OF  PLATINUM  WIRE  RESISTANCE 

the  resistance  value  of  the  platinum  wires  are  obtained  by  measuring  the  voltage 
drop  across  the  desired  wire  and  dividing  it  by  the  current  flow  through  it. 

„  ^PtWire 

^Prme—  I 

‘PifPire 

(1.7233 Volts) 

(0.2831  Amps) 

=  6.0874  Ohms 

1.  Determination  of  Wire  Surface  Temperature 

The  wire's  surface  temperature  are  obtained  from  the  calibration  curves  shown 
in  APPENDIX  A,  wire  1  will  be  shown  as  an  example. 
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"^Ptwirex  ~  53.902  X  Kpiiyifg  —  254.470 
=  53.902  X  6.0874  -  254.470 
=  73.356°C 

2.  Determination  of  Input  Power 

The  input  power  of  the  wires  are  determined  from  the  current  flow  through  the 
wire  times  the  voltage  drop  across  the  wire. 

Powerj  =  ^ ptwirel  ^  ^Ptwirel 

=  I.7233Voltsx0.2831Ohms 
=  0.48787Watts 

3.  Determination  of  Heat  Flux 

The  heat  flux  dissipated  by  the  platinum  wire  is  determined  from  the  inputted 
power  divided  by  the  wire's  surface  area. 


Hflux,  = 


Power 


(0.48787Watts) 


1.5585E-5m" 


=  31,303.44Watts/m^ 


where, 


Aj  —  7r  X  D(  X  T| 

=  TT  X  5E-5  m  X  0.09922m 
=  1.5585E-5'm^ 

C.  DETERMINATION  OF  FLUID  PROPERTIES 
i.  Film  Temperature 

The  film  temperature  is  the  arithmetic  average  of  bulk  fluid  temperature,  and  the 
wire's  surface  temperature. 


T/j/m  — 


Mre) 

2.0 

(54.0"C  +  58.0''C) 


2.0 


=  56.0"C 
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2.  Thermal  Conductivity 

From  Figure  5  of  the  3M  Corporation  Flourinert  Product  Manual  [Ref.  3  :  p. 
15],  the  thermal  conductivity  has  been  determined  to  be: 

(0.6033 -0.001 15 

io. 

=  0.05389Watts/m°C 


'^/llm  —  56  C 

3.  Liquid  Density 

Using  the  expressions  in  Table  4B  and  the  constants  presented  in  Table  4C  of 
the  3M  Corporation  Product  Manual, [Ref.  3:  p.lO],  the  liquid  density  has  been  deter¬ 
mined  to  be: 


Pi  =  (1.740  -0.00261  X  Tyi,J  x  1000 
=  l,593.8kg/m^ 


T/i/m  —  56  C 

4.  Kinematic  Viscosity 

Using  Figure  3  in  the  3M  Corporation  Product  Manual,  [Ref.  3:  p.l3]  and  de¬ 
termining  an  expodintial  curve  fit  valid  from  O^C  to  90'’C  yields: 

V 

A 
B 


=  A  exp( 


B 


T;,,„-f-273.15T 


=  1.203952E  -  8m7s 
=  1, 058.4 109°K 


—  56  C 


V  =  3.000E  -  7m^/s 


5.  Specific  Heat 

From  Figure  4  of  the  3M  Corporation  Flourinert  Product  Manual,  [Ref.  3  : 
p.l4],  the  specific  heat  has  been  determined  to  be: 

Cp  =  (0.241 111  +  3.70374E  -  4  x  TyjJ  x  4186 
=  l,096.IJ/kg’C 


6.  Thermal  Expansion  Coefficient 

Using  the  expression  in  Table  4B  and  the  constants  presented  in  Table  4C  of  the 
3M  Corporation  Product  Manual,  [Ref.  3:  p.lO],  the  thermal  expansion  coefficient  has 
been  determined  to  be: 


_ 0.00261 _ 

^  (1.740  -0.00261  xTyj/J 

=  0.00164  irc 


at  T/„„  =  56°C 

7.  Thermal  Diffusivity 

Is  determined  to  be  at  Ty;,„  =  56°C 


pxCp 

=  3.0848E  -  8m^/s 


8.  Prandlt  Number 

Prandlt  Number  is  determined  to  be  at  =  56°C 


Pr  = 


_v 

a 


=  9.725 
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APPENDIX  C.  UNCERTAINTY  ANALYSIS 
A.  UNCERTAINTY  IN  SURFACE  AREA 

(Od  =  0.003  mm 
=  0.8  mm 
A  =  71 D  L 
where: 

D  =  Wire  Diameter 
L  =  Length  of  platinum  wire 
A  =  Surface  Area  of  platinum  wire 


dA 

3D 

(n  L) 

3A 

dL 

(ttD) 

Therefore, 

co^  = 

[(tt  L  cuof  +  (rt 

II 

3  < 

^  L ;  J 

7  0.003  V  1 

(  0.8  \2~'\ 

_v  0.05  J 

^  99.22  J 

=  0.061  or  6.1% 

A 

03^  =  9.43529E-7m^ 
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B.  UNCERTAINTY  IN  POWER 


(^V2a  =  0.008% +  8mV  [Ref  13,p^.587] 


(^KPtmre  =  0.008% +  8/xV  [Ref  13,/>g.587] 


<y 


«,2a  ±  0.00 1  max  value 


Q  -  ^PiWire  ^  ^ 


PtWire  ^  ^PtlVire 


—  ^PtWire  ^ 


2ft 


2ft 


5Q 


'2ft 


5V 


PtWire 


R 


2ft 


5Q 


'PtWire 


av. 


2ft 


R 


2ft 


dQ 

5R. 


=  -V 


.  +  ■ 


'2ft 


2ft 


Ptiyire''  _2 
‘^2ft 


(Or 


V2ft  \2 

n  ^  (Oy  piiyirg  1  + 

‘^2ft  I 


^Ptmre  \2  .  f  XT  ^^2ft  \2 

^  ^K,2ft  )  +  ( '^P/Wz-e  ^  ^  <^2ft  ) 


R. 


2ft 


1/2 


(Or 


PtWire  \2  ,  (  '^V.2Q.  %  ,  (  «2ft  \l 


V 


+ 


+ 


PtWire 


'2ft 


R 


2ft 


1/2 


The  worst  case  of  uncertainty  will  occur  when  the  voltages  measured  for  Vp,n,.„  and  Vjq 
are  small,  less  than  0.3  Volts. 


(Or 


Q 


(0.008%  +  0.0027%)^  +  (0.008%  +  0.0027%)^  -1  ^  i  9994 100 


1/2 


=  0.0524%  or5.2384E-4 
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C.  UNCERTAINTY  IN  HEAT  FLUX 


Q 

A 


q  =  T- 


_aq_ 

aQ 


5A 


*2 


CO,  = 


(Or, 


1/2 


+ 


CO. 


9  \2 


1/2 


=  [(0.061)  V(5.2384E-4f]'/^ 

=  6.1002E-2or  %  6.1 

D.  UNCERTAINTY  IN  TEMPERATURE 

The  uncertainty  in  the  thermocouple  measurements  (cOtc)  is  determined  by  the  un¬ 
certainty  in  the  platinum  wre  temperature  probe  {(Op,p,^,)  and  the  error  in  the 
thermocouple's  in  matching  the  probe  temperature  (w^t-c)- 

COtc  ~  [(cOpiPrnbef  + 

=  0.400  rc 

E.  UNCERTAINTY  OF  WIRE  SURFACE  TEMPERATURE 

From  the  calibration  formula 

°C 

CO,  =  0.1^ 


=  0.400  rc 
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dTo 


dST 

1 

^^sur/ 

R2a 

^^PlU'ire 

^  Y  n 

^2a 

^T^sur/ 

5R2n 

^Pimre 

®  y 
^2Q 

^Tsur/ 

e>\T 

^2^ 

'^^Ptmre' 

vL 

72 


^Tsurf 


^PtWlre^2a  \2 


<y£  f  +  (I  <^Tof  +  (  ^  "TT^  ‘W^.P/Wre/  +  (  C  (Uja 


•eyptwire )^  +  (i 


Using  data  from  Wire  1  at  the  following  point, 


^  PtWlT, 

=  0.7 1673 Volts 

II 

a 

0.24529  Volts 

^JO  ~ 

1.9994  Ohms  for  wire  1 

^w/r»l  “ 

53.902“C/O 

= 

0.02°C 

This  results  in  an  uncertainty  in  temperature  of, 
^Tsurf,wire\  =  0.854®C 
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APPENDIX  D.  DATA  ACQUISITION  PROGRAMS  USED 


A.  PLATINUM  WIRE  RESISTANCE  CALIBRATION  PROGRAM 

The  following  program  was  used  to  obtain  the  Resistance  value's  for  the  platinum 
wire  specimans  used  in  chapter  2.  This  program  is  written  in  HP  basic  language  for  use 
with  controlling  the  HP-3852A  Data  Acquisition  unit. 


10  DIM  Oht>is(0:61  ) 

20  OUTPUT  709; "INTEGER  I" 
30  OUTPUT  709;"RST  000" 

40  OUTPUT  709:  "US^  000" 

50  OUTPUT  709; "CONE  OHMF" 
S0  OUTPUT  709;"NPLC  1.0" 
70  OUTPUT  709;"NRDGS  G2” 
80  OUTPUT  709; "DELAV  0.1" 
90  OUTPUT  709; "RANGE  10" 


100 

OUTPUT  709; "I  =  99" 

110 

OUTPUT  709; "I  =  1  + 

I" 

120 

OUTPUT  709;''MEAS  OHMF. I" 

130 

ENTER  709;Ohr*i3(»  ) 

140 

Alo=l .E+30 

150 

J=l+J 

160 

Ahi=l .E-30 

170 

A=0. 

180 

FOR  K=l  TO  60 

190 

A=A+0hms( K  ) 

200 

IF  Ohns<K)<Alo  THEN 

210 

Alo=0hois(  K  ) 

220 

ENO  IF 

230 

IF  Ohns(K  )';Ahi  THEN 

240 

Ahi =0hns( K  ) 

250 

END  IF 

260 

NEXT  K 

270 

B=A/60. 

200 

PRINT  "WIRE  SAMPLE  " 

;J 

290 

PRINT  "60  RESISTANCE 

AVERAGE 

IS  ' 

300 

PRINT  "LO  RESISTANCE 

MEASURED 

IS 

’‘;Alo 

310 

PRINT  "HI  RESISTANCE 

MEASURED 

IS 

" ;  Ahi 

320 

IF  J=5  THEN 

330 

GOTO  360 

340 

END  IF 

350 

GOTO  110 

360 

GOTO  370 

370 

DIM  Temps<  0:61) 

380 

K=0. 

390 

OUTPUT  709;"RST  000" 

400 

OUTPUT  709; "USE  000" 

410 

OUTPUT  70S; "INTEGER 

I" 
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420  OUTPUT  709;"CONF  TEMPT" 

450  OUTPUT  709; "NPLC  I .0" 

440  OUTPUT  709;"NRD6S  B2" 

450  OUTPUT  709; "DELAY  0.1" 

460  OUTPUT  709; "RANGE  AUTO" 

470  OUTPUT  709; "1=299" 

480  OUTPUT  709; "1=1+1" 

490  OUTPUT  709;"MEAS  TEMPT,  I" 

500  ENTER  709;Tenp5( *  ) 

510  K=K+1 
520  A=0. 

53G  FOR  L=1  TO  60 
540  A=A+Temps(L) 

550  NEXT  L 
560  Temp=A/60. 

570  PRINT  "THERMOCOUPLE" ;K; "TEMPERATURE" ;TeMp 

580  IF  K=4  THEN 

590  GOTO  620 

600  END  IF 

610  GOTO  480 

620  END 
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B.  MAIN  DATA  ACQUISITION  PROGRAM 

The  following  program  was  used  to  obtain  the  Temperature  value's  for  the  platinum 
wire  specimans  used  in  chapter  3.  This  program  is  written  in  HP  basic  language  for  use 
with  controlling  the  HP-3852A  Data  Acquisition  unit. 


10 
20 
50 
40 
50 
60 
70 
80 
90 
100 
I  10 
120 
150 
140 
150 
160 
170 
180 
130 
200 
210 
220 
250 
240 
250 
260 
270 
230 
290 
500 
510 
520 

j  J0 
540 
550 
560 
570 
580 
390 
400 
410 
420 
450 
440 
450 


DIM  Tenps(0:  16  )  .Tl-(  1  :7  )  ,Ti(  1  ;  I  00  )  .Flu-v-d  :  1  00  )  .Tb(  100  ) 

L=0 

J=0 

L=L+1 

IF  L=81  THEN 
GOTO  1450 
END  IF 

OUTPUT  TOEr-^RST  000" 

OUTPUT  709i"USE  000- 
OUTPUT  709; "INTEGER  I" 

OUTPUT  703;"CONF  TEMPT" 

OUTPUT  709;"NPLC  1.0" 

OUTPUT  709i"NRD6S  17- 
OUTPUT  709 i "DELAY  0.1" 

OUTPUT  709; "RANGE  AUTO" 

OUTPUT  703; "1=239" 

OUTPUT  709i "1=1+1“ 

OUTPUT  703;"MEAS  TEMPT, I" 

ENTER  709iTenps( *  J 

J=J  +  1 

A=0. 

FOR  K=1  TO  15 
A=A+Temps( K  ) 

NEXT  K 
TK< J)=A/15. 

IF  J=7  THEN 
GOTO  500 
END  IF 
GOTO  170 

DIM  Uolts<0: 1  1  )  ,UM  1  : 10) 

J=0. 

Ct=0. 

OUTPUT  709;  "RST  000- 
OUTPUT  709; "USE  000- 
OUTPUT  709; "INTEGER  I" 

OUTPUT  70S;"CONF  DCU" 

OUTPUT  709;"NPLC  1.0" 

OUTPUT  709;"NRDGS  12" 

OUTPUT  709; "DELAY  0.1- 
OUTPUT  709; "RANGE  AUTO- 
OUTPUT  709; "1=99" 

OUTPUT  709; "1=1+1 " 

OUTPUT  709;"MEAS  DCU.I" 

ENTER  709;Uolts<*  ) 

J=J+1 
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4G0  IF  J=5  THEN 
470  GOTO  420 
480  END  IF 
490  Ct=Ct+1 
500  A=0. 

510  FOR  K=1  TO  10 
520  A=A+yolt5(K) 

550  NEXT  K 
540  Vk(Ct)=A/l0. 

550  IF  J=e  THEN 

5B0  GOTO  550 

570  END  IF 

580  GOTO  420 

590  DIM  Uoltw( 0: 1 1  > 

600  J=0. 

610  Ct=5. 

620  OUTPUT  709; "RST  000" 

650  OUTPUT  709; "USE  000" 

640  OUTPUT  709; “INTEGER  I" 

650  OUTPUT  709;"CONF  DCU" 

660  OUTPUT  709;"NPLC  1.0" 

670  OUTPUT  709;"NRDGS  12“ 

680  OUTPUT  709; "DELAY  0.1" 

690  OUTPUT  709; "RANGE  AUTO" 

700  OUTPUT  709; "1=109" 

710  OUTPUT  709; “1=1+1" 

720  OUTPUT  709;“MEAS  OCU.I" 

730  ENTER  709;Uoltw( -► ) 

740  J=J+1 

750  IF  J=S  THEN 
760  GOTO  710 
770  END  IF 
780  Ct=Ct+1 
790  A=0. 

800  FOR  K=1  TO  10 
810  A=A+Uollw(K) 

820  NEXT  K 

850  Uk(Ct  )=A/;0. 

840  IF  J=6  THEN 

850  GOTO  880 
860  END  IF 
870  GOTO  710 

880  Bt  =  (Tk( 1  )+Tk(2  )+Tk(3)+Tk(4>)/4. 

890  Ul=Uk<6) 

900  Il='Jk(l  )/l  .9994 

910  U2='Jk(7) 

920  I2=Uk(2)/2.0076 

950  U3=Uk<8) 

940  15=01.(5)/!  .9997 

950  04=01.(9) 

960  I4=0k(4)/2.0018 

970  05=0k(10) 

980  I5=0k(5)/2.0025 
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990  P1=IHV1 

1000  P2=I2*V2 

1010  P3=I3*U3 

1020  P4=I4*U4 

1050  P5=I5*U5 

1040  R1=U1/I1 

1050  R2=U2/I2 

1050  IR3=U3/I3 

1070  IR4=V4/I4 

1080  iR5=i;5/I5 

1090  Fluxl=Pl/l .5585284E-5 

1100  Flux2=P2/l .5585234E-5  < 

1110  Flux3=P3/l .5710005E-5 

1120  Flux4=P4/l .54605G3E-5 

1130  Flux5=P5/l .5710005E-5 

1140  T1=<53.902123*R1  )-254. 46995 

1150  T2=(53.E66456*R2 >-253.00959 

1160  T3=(53. 141047*R3)-252.729131 

1170  T4=(53.483646*R4 >-252. 848797 

1180  T5=(53.085508*R5 >-260.04883 

1190  Tc=(Tk(5>+Tk(6>>/2. 

1200  PRINT  "  TEMPERATURES" 


1210 

PRINT  " 

AMBIENT  AIR  BULK  FLUID 

CONDENSOR 

SURFACE" 

1220 

PRINT  Tl.,(7>  ,Bt  ,Tc 

1250 

PRINT  " 

"  .L 

1240 

PRINT  " 

WIRE  TEMPERATURE 

CURRENT 

UOLTS  " 

1250 

iPRINT 

<  AMPS  > 

<UOLTS> 

(DEGREE  C>  " 

1260 

PRINT  “ 

1 " .T1  ,11  .U1 

1270 

PRINT  " 

2"  ,T2.I2,U2 

1280 

'PRINT 

"3" .T3,I3.U3 

1290 

'PRINT 

"4''  .T4,I4.U4 

1300 

'PRINT 

"5" .T5,I5,U5 

1310 

PRINT  " 

HEAT  FLUX 

POWER 

RESISTANCE  " 

1320 

'PRINT 

( WATTS > 

(OHMS> 

(WATTS/M2>" 

1 330 

PRINT  " 

1 "  .Fluxl  ,P1  .R1 

1340 

PRINT  " 

2"  .Flux2  .P2  ,R2 

1 550 

'PRINT 

"3"  .FLUX3.P3,R3 

1360 

'PRINT 

"4" ,Flux4,P4,R4 

1570 

'PRINT 

"5"  ,Flux5,P5,R5 

1380 

PRINT  " 

" 

1390 

Tt<L)=T 

1 

1400 

Tb(L>=Bt 

1410 

Flux(L> 

=Flux 1 

1T20  WAIT  100 
1430  GOTO  30 

1440  GOTO  1530  ♦ 

1450  N=1 

1460  PRINT  "  WIRE  TEMP  HEAT  FLUX  BULK  TEMP" 

1470  IF  M=3  THEN 
1480  GOTO  1530 
1490  END  IF 

1500  PRINT  M,Tt(M}  .Flux(M>  ,Tb(M> 

1510  M=M+1 

1520  GOTO  1470 
1530  END 
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C.  NUSSELT  EVALUATION  PROGRAM 

The  following  program  was  used  to  obtain  the  Nusseit  Number  for  the  platinum 
wire  specimans  used  in  chapter  3.  This  program  is  written  in  Fortran  language  for  use 
with  controlling  the  IBM-3803  Main  Frame  Computer. 


FILE:  NUSSELT  FORTRAN  A1 


SUBROUTINE  NUSSELT  (X1,T1,M,Q) 

C 

C  THIS  SUBROUTINE  WILL  CALCULATE  THE  HEAT  TRANSFER  COEFFICIENT  FOR  THE 
C  DATA  OBTAINED  FROM  THE  EXPERIMENTAL  RUNS  BY  USING  THE  CHURCHILL  AND 
C  CHU  CORRELATION 
C 

DIMENSION  X1(30,T1(X),Q(X) 

REALX4  X1,T1,Q,K,NU,NUD 

OPEN ( 29 , FI LE= • RALY • , STATUS= • UNKNOWN ' , F0RM= • FORMATTED • ) 

C 

C  DIAMETER  OF  THE  WIRE  BEING  ANALYSED  IS  NOMINALLY  ,05  MM 
C 

D=5.0E-S 

C 

C  GRAVITY  TERM  OF  9.81  M/S2 
C 

G=9.81 

C 

DO  100  1=1, M 
C 

C  FILM  TEMPERATURE  OF  THE  PHYSICAL  PROPERTIES 
C 

CCC  TF  =  (XKI)  +  TKI)  +  56.0)/2. 

TF  =  (XKI)  +  2.XTl(I))/2. 

C 

C  TEMPERATURE  DIFFERENCE  BETWEEN  THE  WIRE  AND  BULK  FLUID 
C 

CCC  TD  =(X1(I)  -  TKI)  +  56.0) 

TD  =  XKI) 

C 

C  THERMAL  CONDUCTIVITY  OF  FLUID  W/MK 
C 

K  =  (0.6033  -  0.00115XTF)/10. 

C 

C  BULK  MODULUS  OF  THE  FLUID  1/K 
C 

BETA  =  0.00261/(1.740  -  0.00261XTF) 

C 

C  KINEMATIC  VISCOSITY  OF  THE  FLUID  M2/S 
C 

NU  =  1.203952E-8  X  EXP(1058.4109/(TF+273.15)) 

C 

C  SPECIFIC  HEAT  OF  FLUID 
C 

CP  =  (0.24111  +  3.70374E-4  X  TF)X4186.00 
C 

C  FLUID  DENSITY  KG/M3 
C 

RHO  =  (1.740  -  0.00261XTF)X1000.0 
C 


oooo  ooo  ooo  ooo  oo 


THEMAL  DIFFUSIVITY 
ALFA  =  K/CRHOXCP) 

PRANDLT  NUMBER  OF  FLUID 
PR  =  NU/ALFA 

RALYEIGH  NUMBER  OF  FLUID  MUST  BE  10XS(-5<  RA  <  10XX12 

RA  =  (G?(BETA*TDXDXJ(3.)/(NUXALFA) 

PRINT  X,  RA 

NUSSELT  NUMBER  OF  FLUID  NU=HXD/K 

TERM  =  ((0.518x(RAXX(0.25))K(l+(0.559/PR)X*0,6)XX(-5./12.))XX15 
|+(0.1XRA*X(1./3.))XX15.)XX(1./15.) 

NUD=2.0/(LOG(1.+  2. /TERM)) 

I 

THERORICAL  HEAT  FLUX  WATTS/M2 


QCI)  =  KXNUDXTD/D 

^100  WRITE(29,500)  TD,RA,NUD,PR,0(I) 

500  FORMAT  C5(F12.4,2X)) 

CLOSE  (29) 

RETURN 

END 


LIST  OF  REFERENCES 


1.  Bergles  A.E.,  Heat  Transfer  in  Electronic  and  Microelectronic  Equipment, pp  3-60, 
Hemisphere  Publishing  Co.  1990., 

2.  Bergles  A.E.,  and  Kim  C.J.,  "A  Method  to  Reduce  Temperature  Overshoots  in 
Immersion  Cooling  of  Microelectronic  Devices",  paper  presented  at  the  1987 
ASME,  JSME  Thermal  Engineering  Joint  Conference  Honolulu,  Hawaii,  1987. 

3.  3M  Corporation,  Flourinert  Electronic  Liquid  Product  Manual  Commercial  Chemical  Divi¬ 
sion:  St.Paul  Minnesota,  1987. 

4.  You  S.M.,  Simor.  T.W.,  and  Bar-Cohen  A.,  "Experiments  on  Nucleate  Boiling  Heat  Transfer 
with  a  Higlily-Wetting  Di-Electric  Fluid;  Effects  of  Pressure,  Sub-Cooling,  and  Dissolved  Gas 
Content,"  paper  presented  at  the  1987  ASME,  JSME  Thermal  Engineering  Joint  Conference, 
Honolulu,  Hawaii,  1987. 

5.  You  S.M.,  Simon  T.W.,  and  Bar-Cohen  A.,  "Experiments  on  Boiling  Incipience  with  a 
Highly-Wetting  Dielectric  Fluid:  Effects  of  Pressure,  Subcooling  and  Dissolved  Bas  Con¬ 
tent,"  paper  presented  at  the  1987  ASME,  JSME  Thermal  Engineering  Joint  Conference, 
Honolulu,  Hawaii,  1987. 

6.  You  S.M.,  Simon  T.W.,  Bar-Cohen  A.,  and  Tong  W.,  'Experirrental  Investigation  Nucleate 
Boiling  Incipience  with  a  Highly- Wetting  Dielectric  Fluid  (R-113),'  Iniemation  Journal  of 
Heat  and  Mass  Transfer,  v.33,  No.l,pp.l05-117,  1990. 

7.  Marto  P.J.,  and  Lepere  V.J.,  "Pool  Boiling  Heat  Transfer  From  Enhanced  Surfaces  to 
Dielectric  Fluids,"  ASME  Journal  of  Heat  Transfer,  v.l04,  pp.292-299.  May  1982. 

8.  Maddox  D.E.,  and  Mudawar,!.,  "Single  and  Two-Phase  Convective  Heat  Transfer  from 
Smooth  and  Enhanced  Microelectronic  Heat  Sources  in  a  Rectangular  Chaimel,"  ASME 
Journal  of  Heat  Transfer,  v.lll,  pp.  1045-1052. 


81 


9.  Danielson  R.D.,  Tousignant  L.  and  Bar-Cohen  A.,  "Saturated  Pool  Characteristics  of  Com¬ 
mercially  Available  Perfluouinated  Inert  Liquids,"  paper  presented  at  the  1987  ASME,  JSME 
Thermal  Engineering  Joint  Conference,  Honoulul,  Hawaii,  1987. 

10.  Kuehn  T.H.  and  Goldstein  R.J.,  "Correlating  Equations  for  Natural  Convection  Heat 
Transfer  Between  Horizontal  Circular  Cylinders,"  Internation  Journal  of  Heat  and  Mass 
Transfer,  v.l9,  pp.l  127-1 134,  1976. 

11.  Eren  A.S.,  Heat  Transfer  Enhancement  due  to  Bubble  Pumping  in  FC-72  Near  Saturation 
Temperature,  Master's  Thesis,  Naval  Postgraduate  School,  Monterey,  California,  March 
1991. 

12.  Materials  Electronic  Products  Corporation,  MELCOR  FRIGICHIPS  Product  Manual, 
Thermoelectric  Cooler's  Division:  Trenton,  New  Jersey,  1985. 

13.  Hewlett  Packard,  1991  Test  and  Measurement  Catalog,  Pages  586  -  589. 


« 


4 


IT 


82 


INITIAL  DISTRIBUTION  LIST 


1.  Defense  Teclmical  Information  Center 
Cameron  Station 

Alexandria,  VA  22304-6145 

2.  Library,  Code  52 

Naval  Postgraduate  School 
Monterey,  CA  93943-5002 

3.  Mechanical  Engineering  Curricular  OfBce,  Code  ME 
Naval  Postgraduate  School 

Monterey,  California  93940 

4.  Professor  M.D.  Kelleher,  Code  ME/Kk 
Naval  Postgraduate  School 
Monterey,  California  93943-5000 

5.  Professor  Y.  Joshi,  Code  ME/Ji 
Naval  Postgraduate  School 
Monterey,  California  93943-5000 

6.  Naval  Weapons  Support  Center 
Code  6042 

Crane,  Indiana  47522 

7.  LtJg  Ali  Sukru  Eren 
Harzem  Sok  5/3,  80650 
Celiktepe,  Istanbul,  Turkey 

8.  LT  R.A.  Egger 

4830  Rocky  River  Drive 
Cleveland,  Ohio  44135 


No.  Copies 
2 


2 


1 


2 


1 


1 


1 


2 


* 


* 


83 


